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We study the Josephson effect between weakly coupled d-wave superconductors within the qua- 
siclassical theory, in particular, the influence of interface roughness on the current-phase relation 
and the critical current of mirror junctions and 45° asymmetric junctions. For mirror junctions 
the temperature dependence of the critical current is non-monotonic in the limit of low roughness, 
but monotonic for very rough interfaces. For 45° asymmetric junctions with a linear dimension 
much larger than the superconducting coherence length we hnd a sin(2i£>)-like current-phase rela- 
tion, whereas for contacts on the scale of the coherence length or smaller the usual sin if -like behavior 
is observed. Our results compare well with recent experimental observations. 

PACS numbers: 74.76.Bz,74.50.+r 



I. INTRODUCTION 

It is commonly believed that high-T c superconduc- 
tors are described by an order parameter which domi- 
nantly exhibits d-wavc symmetry This was especially 
established by SQUID-like experiments, e.g. the corner- 
SQUID experiments by Wollmann et ali or the tri- 
crystal experiment by Tsuei et al^, which provide a 
phase-sensitive test of the order parameter. This ob- 
servation was followed by numerous theoretical studies 
of the Josephson effect between d-wave superconductors 
with rather unexpected results; a review can be found in 
ReisMZ. 

For instance a mirror junction, where the order pa- 
rameter on the left and right hand side is rotated by 
the same angle but in opposite directions, should ex- 
hibit a non-monotonic temperature dependence of the 
critical current2&i; this behavior is related to a tran- 
sition to a 7r-junction state at low temperatures. For 
the 45° asymmetric junction, in which the order pa- 
rameter is rotated only on one side by 45°, a sin(2<y9)- 
like current-phase relation was predicted^, as the lead- 
ing contribution to the tunnel current vanishes due to 
symmetry. Experimentally the temperature dependence 
of the critical curren1i2i2iiS and the current-phase rela- 
tionii^u3 f artificially produced grain boundaries with 
a well defined lattice orientation were determined. For 
mirror junctions the predicted non-monotonic behavior 
of the temperature-dependent critical current was found 
by Il'ichev et alii^, whereas in other cases a monotonic 
behavior as in usual s-wave superconductors was re- 
ported&Siifi. In the experiment on an asymmetric junc- 
tion^ not all samples showed the predicted sin(2</?)-like 
current-phase relation. 

As the d-wave order parameter is sensitive to disorder 
scattering it is believed that the origin for these differing 
results is interface roughness: One possibility is the ex- 
istence of facets with spatially varying orientation which 
occur at the interfaces on a /xm-scal o 14 ' 15 . This kind of 



disorder can theoretically be described by a model con- 
sisting of randomly oriented mirrors 1 ^!. On the other 
hand roughness on the atomic scale is present as well. 
In theoretical studies this is a non-trivial problem which 
we will tackle in this work. In previous studies inter- 
face roughness was modeled by a thin dirty layer next 
to a specular interface^* 1 ^*^. We use an alternative ap- 
proach which allows us to study, at least in principle, also 
individual realizations of the disorder; this might be of 
importance for mesoscopicly small junctions. 

For the calculations we will use the quasiclassical the- 
ory of superconductivity22i2i which is valid for spatial 
variations on scales which are large compared to the 
Fermi wave length, A^; this is marginally fulfilled in the 
Cu02-planes of high-T c materials as £o/ ~ 5 (£o : su- 
perconducting coherence length at zero temperature in 
the Cu02-planes). Near interfaces the quasiclassical the- 
ory must be supplemented by boundary conditions. A 
specular interface can be described by Zaitsev's^ bound- 
ary conditions. We use a more general scheme as sug- 
gested by Shelankov and Ozana2£; this approach is ex- 
tended in order to find suitable boundary conditions for 
interfaces. The properties of the interface, e.g. roughness 
or transparency, are incorporated by a scattering ma- 
trix. To describe an interface without regular structure 
we phenomenologically choose random scattering matri- 
ces. A similar approach has already been successfully 
applied to examine tunnel junctions between a normal 
metal and a d-wave superconducto»2iiS£. 

The paper is arranged as follows: In the next section 
we briefly introduce the quasiclassical theory. We then 
describe in detail the boundary conditions used in this 
work, and discuss our choice of the scattering matrix and 
its properties. In section ITTT1 we first present our results 
for mirror junction: We study the temperature depen- 
dent critical current for interfaces with varying rough- 
ness and make a quantitative comparison with experi- 
mental data. We furthermore examine the transition to 
a 7r-junction, where a spontaneous current parallel to the 
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junction is present. Secondly we present the results for 
45° asymmetric junctions, where additionally to the av- 
erage quantities also statistical fluctuation of the critical 
current are examined. In both cases we discuss in detail 
the influence of interface roughness, which might lead to 
clear modifications of the results for the clean case. Con- 
cluding remarks are given in section llVl 



II. METHOD 

In this section we introduce the main ingredients for 
our phenomenological model of rough superconducting 
contacts: The superconducting state is described by the 
quasiclassical theory, which is supplemented by boundary 
conditions to take into account interface scattering. We 
first recall the basic formalism. 



A. Theory of quasiclassical Green's functions 

The quasiclassical matrix Green's function is deter- 
mined by the Eilenberger equation which in thermal equi- 
librium reads 



t 3 E + ef 3 v F ■ A(r) + iA(p F , r), g(E, p F ;r) 
+iv F ■ d r g(E, p F ]r) = 



(1) 



Here % are the Pauli matrices in Nambu space. The 
physically relevant solution is fixed by the normalization 
condition 

[g(E,p F ;r)} 2 = i. (2) 
The (spin-singlet) order parameter has the matrix form 



A(p F ,r) 







A( P F,r; 



v A*( PF ,r) 
and must be determined self-consistently via 



(3) 



A(pf,r) = - 7 rA/'oT ]T (V(p F , p' F )g(iE n , p' F ; r))^ , 

£„|<_E C 

(4) 

where E n — nT(2n + 1) are the Matsubara energies 
(h = he = 1) and Wo is the density of states (DOS) 
at the Fermi level in the normal state. For convenience 
we assume a cylindrical Fermi surface; (. . .) denotes 
the Fermi surface average. In the BCS approximation 
the pairing interaction is given by its strength V, its cut- 
off energy E c , and its direction dependence which can be 
expanded for a d-wave order parameter as follows 



V(pf,p' f ) = Vt](pf)v(p'f) 



with 



v(Pf) = (Pf.x -Pf v )/p 



F ■ 



(5) 



(6) 



For a homogeneous system and temperature T = one 
finds A(p F ,r) = hA r)(p F ), with A « 2.14T C . 

Having solved these equations for the Green's function, 
various observables can be calculated. In this work the 
current density (in thermal equilibrium) is of particular 
interest: 



j(r) 



■A/" T Tr[f 3 (v F ff(iE n ,p F ;r)) 



(7) 

Another important quantity is the DOS which is given 
by 

M(E,r) = (M(E,pF,r)) pp , (8) 
with the angle-resolved DOS 

M(E, p Fl r) = ^M Rc{Tr [r 3 g(E + iO + ,p F ;r)]} . (9) 

In order to solve the Eilenberger equation it is useful 
to parameterize the Green's function as suggested in 
RefsM-21 



f 
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1 - ah 



1 + ab -2a 
2b -(1 + ab) 



(10) 



Note that the normalization condition is fulfilled by con- 
struction. The ansatz (|1U|> yields the Riccati type equa- 
tions for the amplitudes a(E, p F ;r) and b(E, p F ;r): 

v F ■ d r a = A*a 2 + 2i(E + ev F • A)a - A , (11) 
v F ■ d r b = Ab 2 -2i(E + ev F ■ A)b- A* . (12) 

The (numerical) evaluation of these amplitudes along 
classical trajectories r(A) = ro + \v F /v F is more con- 
venient than solving the Eilenberger Eq. together 
with the normalization condition Eq. In the case 

Im[£T] > 0, i.e. for the retarded Green's function or the 
Matsubara Green's function with E n > 0, the integra- 
tion of Eq. (|llf> is stable in positive v F -direction, and 
Eq. i|12|) is stable in negative v F -direction; the directions 
must be reversed for Im[_E] < 0. In the cases considered 
here the trajectories start and end in the bulk (A = ±oo) 
where the order parameter has no spatial dependence; 
then Eqs. I|llf) and H12|) yield the homogeneous solutions 
(d r a = 0, d r b = 0) 



a(£;,p F ; A 



b(E,p F ;X 



-co 



-iA_(p F ) 



£ + v /£ 2 -|A-( PF )| 2 ' 

iA* + (p F ) 
E+^E2-\A + (pf)\ 2 ' 



(13) 



(14) 



where A±(p F ) is the order parameter for A — > ±oo. 
Starting from these initial values Eqs. I|ll|) and (|12|) must 
be integrated in positive and negative directions, respec- 
tively. It is then possible to construct the Green's func- 
tion on the trajectory via Eq. (|TU|l . 

In the following subsection we show that this param- 
eterization is a good starting point for implementing 
boundary conditions which are necessary to study inter- 
faces. 
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B. Boundary conditions 

The amplitudes a with pointing towards the in- 
terface (in-direction) can be calculated by integrating 
Eq. with the initial value in the bulk given by 

Eq. l(T3"|) ; analogously the 6's with pointing away 
from the interface (out-direction) can be calculated via 
Eqs. 1(12(1 and 1)14(1 . The boundary conditions provide the 
a's in out-direction and the o's in in-direction directly at 
the interface; starting with these values the integration 
of Eqs. 1(11(1 and 1(12(1 . respectively, yields the a's in the 
out- and the 6's in the in-direction. 

For simplicity we consider only a finite number n of 
trajectories on each side of the junction, i.e. P' s Fin / out — * 

PFin/out = 1' ■ • ■ ' n ' s = '/ r: left /right); we choose the 
numbering of the directions such that (see Fig. ^) 



l,i 
PFin,y 



PFout, y 



PFb 



r.i 
PFo 



(15) 

We denote the amplitude a on the i th in/out-trajectory 
on the side s of the junction by af n / out , an< ^ ^in/out ana l~ 
ogously. To formulate the boundary conditions we follow 
RefiSS by introducing the determinant T> as follows: 



V 



= det (l - SaS+b) ; (16) 
the 2n x 2n diagonal matrices a and b read 

a = diag {a l ^},{a r ^} , (17) 

b = dfcg (18) 

where the a's and 6's are the amplitudes at the interface 
(x = 0). The unitary 2n x 2n scattering matrix S, which 
contains the physical properties of the interface, has the 
block structure 

2rl err (19) 

where the elements are n x n matrices. As described in 
detail in2£ we define the quantities 

dV 

af ' =0 



V 



b"&=0 



da in 

dV 



b s A=o 



(20) 
(21) 









y ~x 







FIG. 1: At an specular interface the Green's functions on the 
four trajectories with the same parallel momentum pf, v are 
coupled coherently. 



The boundary conditions yield the unknown values of a 
and b at the interface via 



B s { 



(22) 



To apply these boundary conditions to interfaces, we 
must ensure current conservation across the junction. As 
discussed in Refi2& the above boundary conditions yield 
the following conservation law: 



4 = 1 



n 

=£[mp r /out)-mp£D] 



(23) 



where the contribution of each direction is given by 

~f(E,p F ) =Tr[T 3 g s (E,p F ;Q)}. (24) 

On the other hand, we impose current conservation per- 
pendicular to the junction (see Eq. J2J) 

(v F , x [j l (E,p F J-j l (E,p l Fout )] ) 



PiTu 



(v F , x [f(E,p Fia )-~f(E,p Fout )]) 



(25) 



PFin 



Note that the directions p^j n and p s Font are related by 
Eq. (fL5)l . see Fig.^ To ensure current conservation per- 
pendicular to the junction, we construct the grid of the 
discrete directions so that the term v Fx = v F p FtX /p F is 
already taken into account. This is guaranteed by choos- 
ing the grid as 



( cos di 
/ - COS 

PFout = PF 



■ COS 1?, 
PFin=PF\ sin1? . 



sinz9i = 



2i 
n+l 



sini? 

-l, ; = :i. 



r i I COS 

PFout = PF 



sint?i 



(26) 



In other words this grid takes into account that the rate of 
scattering events at the interface decreases with smaller 
angles of incidence. 

The probability of scattering from the direction p^; n 



into an interval [p^jP^ 1 ] (s, s' = l/r) reads 



Ps.tfj-ttyAtii = \S- 



s' s\2 I 



(27) 



with Ai?j = -di — Using the non-equidistant grid 

as defined in Eq. 1(26(1 for a large number of scattering 
channels, n 3> 1, we find 



P s < s (^ -COS^ s | 2 



(28) 



In the next step one has to find the appropriate scat- 
tering matrices for a given physical realization. In the 
following subsection we study the case of a specular and 
an irregular rough interface. 
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C. Scattering matrix 

The most simple case is a specular interface where only 
trajectories with the same parallel momentum are cou- 
pled. Then the block matrices that form the scattering 
matrix in Eq. (|19fl have diagonal form: 



S" = -S rr = R = diag a/1 -T(tfi) 

S lr _ S lr = r 



diag v / T^F l 



(29) 



in particular, for a 5-like boundary potential the 
direction-dependent transparency reads 2 ^ 



T cos 2 d 
I -T sin 2 



To e [0,1]. (30) 



For a normal metal this leads to the resistance of the 
tunnel-barrier 



R 



_! _ 4:Ae 2 Af VF 



N 



3tt 



To, 



(31) 



where A is the area of the contact. 

In the superconducting case the boundary conditions 
for a specular interface can be solved for each direction 
individually, with the resuhSiSi 

a s,i = a in a in 6 out ) + a in 7^(1 - a in 6 out ) 

^(1 - a in & out) + - a i„ »out) 

,sA o out T,(l - a in 6 out ) + b out Ki{l - a in b out ) 



T i( l ~ a in & out) + - afn^out) 



(33) 



where s = Z/r, s = r/l, respectively, and T = T(d{), 
1*4 = 1-%. 

A rough surface, on the other hand, is described by a 
scattering matrix that does not only depend on the trans- 
parency but also on a roughness parameter. Our starting 
point to construct such a unitary scattering matrix is 31 



S = 





) 


I 




(U 3 


) 


o 1 













(34) 



with R and T as defined in Eq. I|29[l . which contain the 
information on the transparency. The unitary matrices 
Ufe describe the interface roughness; for = 1 the 
scattering matrix for a specular interface obviously is re- 
covered. 

We will focus on interfaces without regular structure. 
To take into account the statistical character of rough 
interfaces we choose the Ufc to be unitary random ma- 
trices 2 ^, 



U fe = cxp{iH fe }, 



(35) 



where each is a random matrix with Gaussian corre- 
lations: 



T 

{Hk,ij) — , {HkAj*Hk'a>j') — -^Su'Sjj/Skk' ■ (36) 



s 




FIG. 2: The weight for specular scattering, \u\ 2 , decreases 
with growing r; for r > 3 specular scattering is almost com- 
pletely suppressed (|it| 2 <C 1). 



Here (...) denotes the disorder average. In the following 
we use the abbreviations 



|«(t)| = (\UkM 2 ) 



HT)\=(\U k ^\ 



(37) 



to describe the average scattering probability; due to the 
unitarity it follows that |u| + (n— l)|u| = 1. Obviously for 
t = we find |it(r)| = 1, whereas for increasing (positive) 
r, we obtain a reduced value, |w(r)| < 1, as can be seen 
in Fig. |21 A detailed discussion can be found 

Using Eq. (|28J) we find the average transmission prob- 
ability density for scattering from the direction p^; n on 
one side of the junction to p^ ut on the other side, 



cos 



{nM 2 T(^)c^+«(i-M) 2 



(P ss {d 3 - di)) 

+H(i-H)[T( 1 ? i )+m-)](i-%)} . 



with 



for To < 1 
for T = 1 



(38) 



(39) 



Furthermore the reflection probabilities (P ss ) can easily 
be obtained by substituting T — > (1 — T) and n — > (1 — fs) 
in Eq. (|3"%)l . The continuum limit of Eq. (|3"%|l reads 

(P aS (& -> #)) = \u\ 2 T(d')5{d' - 0)+ 

+ - m) 2 + i«i(i - m)fw + n<&)]} ; 

(40) 

note that the weight of specular scattering is determined 
by the quantity \u\ 2 : For a specular interface (r = 0) it 
is u.| 2 = 1, whereas in the very rough limit (r 3> 1) we 
find \u\ 2 (see Fig.EJ. 

So far we discussed the boundary conditions at a single 
point of an interface for a discrete number of directions 
pp. A real interface will have a finite cross section and 
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the directions will be continuous. Therefore the question 
arises how these boundary conditions can be applied to 
real interfaces. 

First we consider the discretization of the Fermi sur- 
face. A finite number n sets the typical angle, $ c « ir/n, 
over which the scattering probability can vary. For small 
?9 C the fluctuations of the scattering probability with 
varying directions are averaged out more effectively by 
the Fermi surface average that must be taken to evalu- 
ate observables such as the current or the order parame- 
ter. Therefore the self-averaging character of the physical 
quantities increases with decreasing d c ; i.e. the statisti- 
cal fluctuations of physical quantities are diminished for a 
small d c . It is important to note that only the statistical 
fluctuations of physical quantities depend on the value of 
i? c , but that the mean values are independent of i? c , as 
long as d c <C tt. Throughout the following calculations 
we choose n = 40. 

Second we comment on the number of different real- 
izations of the scattering matrix that must be taken into 
account. In the original formulation of boundary con- 
ditions by Shclankov and Ozana^ the spatial resolution 
is limited by the Fermi wave length, 1/pf', therefore the 
number of different realizations should be of the order 
App, where A is the area of the contact. On the other 
hand, we expect that quantities like the current density 
or the order parameter, which are of interest to us, vary 
on the scale of the coherence length, £o- Therefore we 
choose one effective scattering matrix S to describe an 
area £q °f the whole contact. As the order parameter 
relaxes on the coherence length it is reasonable to treat 
each of these areas individually neglecting the influence 
from other part of the contact; this approximation con- 
siderably simplifies the self-consistent treatment of the 
quasiclassical theory. Finally the current perpendicular 
to a junction of area A is calculated via 

I x = / dA^(r)«£oEj^' (41) 

A 1=1 

where each is evaluated for one particular realization 
of the scattering matrix. For large junctions (A ^> £q) 
the current is given by the average over many scattering 
matrices; for smaller junctions (A of the order of £g or 
even smaller) statistical fluctuations become relevant. 

III. RESULTS AND DISCUSSION 

We consider two kinds of Josephson junctions. First 
we discuss so-called mirror junctions where the order pa- 
rameter on both sides is rotated by the same angle a but 
in opposite directions (see Fig. UJ. Second we examine 
the 45° asymmetric junction where the order parameter 
is rotated only on the right hand side (see Fig. . 

We recall the general expansion of the current-phase 
relation between two superconductors, 

I x (T, ip) = h (T) sin ip + I 2 (T) sin(2^) + . . . (42) 



with the superconducting phase difference, ip, between 
the left and the right hand side. It is also worth men- 
tioning that the current-phase can be obtained from the 
free energy of the junction via 

4CZ» = 2|e|^CZ». (43) 

Besides their temperature dependence, the contributions 
I k scale as 

h oc % k (44) 

for small transparencies. For s-wave superconductors 
and small enough Tq only the linear term in the trans- 
parency, I\ , is relevant. In the following we will show that 
this needs not be the case for rf-wave superconductors. 

A. Mirror junction 

In this section we will consider mirror junctions, which 
are schematically illustrated in Fig. [3] For a junction 
with purely specular scattering, a particular symmetry 
is present, which reads 

A l (p l Fin ) = A r (p r Fin ), 

A i (pU) = A(pU. 

This means that for a given incoming trajectory the or- 
der parameter is identical for the transmitted and the 
reflected quasi-particle. All directions can therefore be 
divided into two classes: (i) For some directions the 
sign of the order parameter is the same for all involved 
trajectories (A'(p^ in )A i (p^ out ) > 1). The behavior of 
these directions is similar to that of usual s-wave su- 
perconductors: Their contribution to I\ is positive and 
finite for T < T c . (ii) For the other directions the 
sign for the in- compared to the out-trajectories changes 
(A'(p^ in )A ; (p^ out ) < 0). This has two crucial conse- 
quences: These directions have a negative contribution 
to I\ which moreover diverges in the tunnel limit as 1 /T 
for T — * 0. The reason for this divergence is the existence 



SC I sc 




FIG. 3: For mirror junctions the order parameter on the left 
and the right hand side of the junction are rotated by the same 
angle a but in opposite directions. For specular scattering the 
order parameters on both in-trajectories and out-trajectories, 
respectively, are equal. 
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of Andreev bound states near zero energy at the inter- 
face as the reflected quasi-particles acquire a sign change 
of the order parametenS^; the divergence is cut-off by a 
finite transparency, which leads to a shift of the bound 
state to finite energies, and by disorder, which leads to a 
broadening of the level. 

Altogether the s-wave-like contributions dominate for 
high temperatures, which, for a fixed phase difference 
< tp < 7r, leads to a positive current, whereas for low 
temperature the anomalous contributions are enhanced, 
which results in a negative current. In other words, 
the ground state of the junction shifts from tp = In- 
junction) to tp — 7r (7r-junction) when decreasing the 
temperature. This transition occurs at the temperature 
T ff , where both contributions cancel; at this tempera- 
ture, the leading sin ip-Uke contribution vanishes (I\ = 0) 
and the junction is dominated by higher order terms. As 
the amount of directions preferring a 7r-junction increases 
with a growing angle a, also the temperature T n increases 
with a. 

In the following we will discuss the influence of in- 
terface roughness on the temperature dependence of the 
critical current. We calculate the current contributions 
Ii and I2 as defined in Eq. I|42(l : higher order contribu- 
tions, Ik, k > 2, can be neglected as we consider only 
small transparencies. We determine the critical current, 
J c , from the absolute maximum of the current phase rela- 
tion. In the graphs a 7r-junction behavior is indicated by 
a negative value. It should be mentioned that in exper- 
iments the absolute value, |/ c | is measured, thus a min- 
imum in the temperature dependent critical current is 
observed at T = T^. We focus on the case a = 22.5°; for 
the presented results the transparency is set to Tq = 0.01. 
The average current is evaluated by using 20 realizations 
of the scattering matrix; the statistical error is less than 
10%. 

First, we concentrate on the clean case. There, as men- 
tioned above, the directions preferring a 7r-junction and 
those preferring a 0-junction compete: From the tem- 
perature dependence of the critical current, as shown in 
Fig. [3 we find a 7r-junction behavior below T v w 0.36 T c , 
whereas for high temperatures the 0-junction state is fa- 
vored. Recalling Eq. (|43|1 . the shift of the ground state 
from tp — to tp — 7r can be seen in the current-phase re- 
lation for decreasing temperature near T v (see Fig. 01 . 
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FIG. 5: Critical current for a mirror junction with 7o = 0.01 
and a = 22.5°. Negative values indicate a 7r-junction state. 
Near T = T n the plotted curves make a jump from a positive 
to a negative value and hence the critical current stays finite. 
Due to the small transparency the jump is not visible in this 
plot but in Fig. |S] 
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FIG. 4: For a 45° asymmetric junction the order parameter 
is tilted on one side only. 



FIG. 6: Current phase relation of a mirror junction without 
roughness (To — 0.01, a = 22.5°) near the temperature where 
it becomes a 7r-junction (T = T 7I ±k- 10~ 3 T C , k = 0, 1, 2, 3, 
increasing T from bottom to top). The asterisks mark the 
points which define the critical current. Note that the criti- 
cal current remains finite at T = T w due to the sin(2<y9)-like 
contribution. 



In particular a finite critical current can be observed 
even for T = T w , since the contribution I2 remains fi- 
nite, whereas I\ vanishes. In Fig.[3]the jump in the cur- 
rent &tT = T 7T is too small to be visible; the reason is the 
tiny transparency (Tq = 0.01) we used in this calculation. 
Moreover it can be seen in Fig. that the current-phase 
relation at is sin(2<^)-like, which is related to the de- 
generacy of the ground state. The ground state exhibits 
a non-trivial phase difference tp = tt/2, 37t/2 which leads 
to a current parallel to the interface; notice that also the 
anisotropy of the order parameter plays a crucial role 
such that the contributions with opposite y-component 
do not cancel each other. The parallel current is pre- 
sented in Fig. where we use jo = 2eNoVF&-o as a unit 
for the current density. We should mention that the cur- 
rent scales with the transparency of the contact, T$. For 
other orientations of the order parameter, a < 22.5°, the 
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FIG. 7: Ground state current in y-direction at a mirror junc- 
tion with a = 22.5° and % = 0.01 at T = T w , i.e. ip « tt/2. 
The current density is increasing with the roughness of the 
interface. 
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FIG. 8: Critical current of mirror junction with different ori- 
entation for very rough interfaces (r = 4). The temperature 
dependence of I C (T) /I C (0.05T C ) is monotonic (within the nu- 
merical error) and its angle dependence is only weak. 



same qualitative behavior should be observed but with a 
reduced T w . 

Obviously, roughness suppresses the 7r-junction behav- 
ior (see Fig. EJ. With increasing r, the temperature 
decreases until the transition disappears for very rough 
interfaces. This can be understood as follows: The nega- 
tive current contributions are carried by Andreev bound 
states with E < C 5 ; as they are broadened by disorder 25 
their negative current contribution is partially canceled 
by those bound states with E > 0. Due to this suppres- 
sion of the 7r-junction contribution at rough interfaces, 
the normal 0-junction contribution becomes dominant 
also for lower temperatures; this contribution is more sta- 
ble against roughness as it is carried by bound states at 
E > — |A|£. For large enough disorder we finally find 
I\ > in the whole temperature range. When further 
increasing the roughness the critical current is enhanced 
until the negative contributions vanish completely Then, 
the temperature dependence of the critical current, i.e. 
/ C (T)// C (0.05T C ), has a monotonic shape almost as in 
the s-wave case. For a very rough interface the quantity 



I C (T) / I c (Q.0bT c ) depends only weakly on the orientation, 
a (see Fig. EJ). 

Also with disorder the ground state at T = T v (as long 
as T-k > 0) exhibits a finite phase difference, which leads 
to a spontaneous current parallel to the interface. As 
shown in Fig.0the current density increases rapidly with 
growing interface roughness. This can be understood as 
follows: In the clean case the directions preferring a 0- 
junction and those preferring a 7r-junction contribute to 
the current with opposite sign; this leads to a partial 
cancellation also of the parallel current. As discussed be- 
fore, roughness suppresses the 7r-junction behavior; as a 
consequence the cancellation of the current contributions 
is less effective, hence the current increases with rough- 
ness. It should be mentioned that for this type of junc- 
tion (a = 22.5°) the modification of the order parameter 
is small and plays only a minor role. 

The suppression of the 7r-junction behavior by disor- 
der has been reported earlier for alternative models of 
interface roughnes a 6 ' 19 . 

In the following, we will compare our results for 
the critical current with the experimental data of 
r^ e f S| 8 4 9 4 i4 i 32_ There a monotonic temperature depen- 
dence of the ratio I(T)//(0.05T C ) has been reported, in 
agreement with our observations for very rough inter- 
faces. The quantitative agreement of these results with 
the data reported inS for YBCO-junctions is reasonable 
for small angles, a, whereas for larger angles it is quite 
poor; this can be seen when comparing the values of RnIc 
at T = 4.2 K w 0.05 T c (see table 0. The discrepancy 
between our results and the experimental data might be 
related to facets on /xm-scale at the interfac o 15 ! 34 , which 
have an increasing influence for larger angles a, as more 
facets yield a negative current contribution 34 . An aver- 
age over all facet orientations then leads to a diminished 
total current. 

Recently, for a small junction (width 0.5/mi) with 
a = 22.5° a non-monotonic temperature dependence of 
the critical current has been reported 1 ^ which, as already 
discussed, appears due to the transition to a 7r-junction 
for low temperatures. The minimum of the critical cur- 
rent was found for = 12 K = 0.13 T c . The trans- 
parency is quite high, as can be seen from the large 
contribution. The measurement is in qualitative agree- 



12° 
18° 
22.5° 



AR N (flcm 2 
Refj 

5.4 x 10" a 

1.5 x 10" 8 
1.2 x 10" 8 



,-j-cxp 



i?jv/ c (mV) 
Ref£ 



R N I c (mV) 
calculated 



0.024 1.3 
0.008 0.75 
0.011 0.13 



3.6 
2.6 
1.9 



TABLE I: Comparison of experimental results for RnIc 
with calculated values. The transparency in the experi- 
ment, Tq* p , is estimated via Eq. I|31|l leading to the val- 
ues ARn given in the table, and vf = 4.34 x 10 6 cm/s, 
./Vo = 2.13 x 10 22 /cm 3 eV^. For the calculations we used 
T = 0.01 (tunnel- limit), T = 0.05 T c , and t = 4. 
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FIG. 9: Average current I2 for an asymmetric junction with 
Tq — 0.01 and varying strength of the roughness; interface 
roughness suppresses the current 1?.. 
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FIG. 10: Current-phase relation for single realizations of the 
disorder for an asymmetric junction with T = 0.1T C . The 
transparency is 7o = 0.01 and the roughness is given by the 
same parameters r as in Fig.|5] 



ment with results obtained from our model for Tq > 0.2 
and t > 0.4. It is also worth mentioning that this tran- 
sition was observed for some of the samples only, which 
might be due to the sensitivity of mesoscopic junctions 
to fluctuations of the interface properties. Note that here 
facets are of minor importance, as the width of the junc- 
tion is comparable to the typical facet size; so it is pos- 
sible to have junctions with well-defined orientation. 

In summary, we showed that the behavior of mirror 
junctions can be modified drastically by interface rough- 
ness: If the roughness is weak and the angle a is large 
enough, a transition to a 7r-junction at low temperatures 
can be observed; near the transition temperature, T„, 
the junction has a sin(2<y9)-like current-phase relation. 
For larger roughness the 7r-junction behavior is destroyed 
and the quantity I C (T) / I C (0.05T C ) depends only weakly 
on the orientation a. When comparing with experimen- 
tal data, we find qualitative agreement with our theory. 
However, especially the influence of large facets has to 
be taken into account to determine the absolute value of 
the critical current correctly. 



B. 45° asymmetric junction 

For a clean 45° asymmetric junction, the particular 
geometric symmetry is responsible for a vanishing leading 
current contribution I\. The symmetry operation y — > 
—y leads to a phase-shift ip — > ip + 7r, while on the other 
hand, the current in rc-direction must be invariant. We 
therefore find the following condition for the current: 

I*(<p) =I x (<P + it) & hk-i = 0, keN, (46) 

and the current-phase relation takes the form 

I x if) = h sin(2^) + 7 4 sin(4^) + . . . . (47) 

This means that in the absence of disorder the current is 
dominated by the contribution /2sin(2y;). Asymmetric 
junctions therefore have two degenerate ground states 
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FIG. 11: Standard deviation of the current I\ for an asym- 
metric junction with To = 0.01 and varying strength of the 
roughness, i9 c ~ 7r/40; this current contribution is increasing 
with the roughness. 



at ip — 7r/2,37r/2. For this reason such junctions are 
discussed as possible realizations of quantum bits^Si. 

Considering rough interfaces we will first examine the 
average values of 1\ and I2; later we will also discuss 
the statistical fluctuations and their relevance for physi- 
cal realizations. We present data for a tunnel-junction 
(Tq = 0.01). Moreover in our approach the statisti- 
cal fluctuations depend on the angle i? c which is set to 
■d c ks ir /40 (i.e. n — 40 directions were taken into ac- 
count). 

Since, in our model for rough interfaces, the symmetry 
y — > — y is still present on average, it follows that (Ix) — 
0; i.e. on average the tunnel-current exhibits a sin(2<ys)- 
like current-phase relation. 

The temperature dependence of (I2) is shown in Fig. El 
for various roughness parameters. In the clean case for 
T — * and T — > 0, we find I 2 oc 1/T due to the zero 
energy bound state at the interface; this state occurs as 
a result of the sign change of the order parameter for re- 
flected quasi-particles on the right hand side. But finite 
roughness also leads to a suppression of this contribution 
and the average value (I2) decreases with growing inter- 
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FIG. 12: Average ground state (i.e. ~ tt/2) current density 
in y-direction at an asymmetric junction with To = 0.01. 



face roughness. As (Ii) = the ground state of the junc- 
tion has a non-trivial phase difference <p = 7r/2, 37r/2, and 
hence the time-reversal symmetry is broken. Thus a cur- 
rent parallel to the interface exists (see Fig. I12|) whereas 
no current in rr-direction is present; as in the previous 
case the parallel current scales with transparency. The 
parallel current is carried by bound states at E < 0, 
which exist on both sides of the interface due to the fi- 
nite transparency^; the bound states with E > would 
contribute to the current in opposite direction, but they 
are not occupied. Interface roughness generally leads to a 
suppression of the parallel current since the bound states 
are broadened, and their current contributions tend to 
cancel each other. Only in a narrow region on the right 
hand side of the interface, of the order of the coherence 
length, the parallel current is enhanced, which is related 
to the considerable increase of the tilted order parameter 
at a rough interface 2 ^. 

The average critical current of 45° asymmetric junc- 
tions was considered before for the clean&SLS as well 
as for the rough case using the thin dirty layer model^. 
The results are similar to ours. The spontaneous parallel 
current was studied for the clean case2&. Roughness was 
also considered for a completely transparent interface 40 , 
where a reduction of the parallel current was found for 
increasing roughness. The sin(2<£>)-like behavior has been 
observed in an experiment with a YBCO-junction by 
Il'ichev et ali£. 

So far we assumed that the junction is well described 
by the average values of the current, but when consid- 
ering a particular realization of the disorder (i.e. of the 
scattering matrix) the symmetry stated in Eq. 146|) is 
no longer valid, and a finite contribution 1\ is expected. 
To confirm this assertion we consider the current-phase 
relation for one particular scattering matrix at the tem- 
perature T — 0.1 T c ; the result is shown in Fig. ^| The 
sin(2</j)-behavior, observable in the clean case, is clearly 
modified by the roughness and a strong sin tp-part addi- 
tionally occurs; for r > 0.4 the current-phase relation is 
already dominated by this contribution. In order to study 
the statistical properties of the current more systemat- 



ically we consider the standard deviation of the contri- 
bution I\, namely ((A/i) 2 ) 1 / 2 ; its temperature depen- 
dence is shown in Fig. ^2 Even for small roughness the 
quantity ((A/i) 2 ) 1 / 2 may already be of the same order of 
magnitude as {12), in particular for higher temperatures 
(compare Figs. l9l and II 1|) . The reason is the different 
dependence on the transparency of both contributions, 
I\ oc 7o and 1% oc 7^ 2 , so that for T <c 1 the contribution 
l\ dominates as soon as it is allowed by symmetry. 

Recalling the physical meaning of the averaging pro- 
cess (see Sec. Ill C|l we conclude that the statistical fluc- 
tuations should be relevant for small junctions (A of 
the order of £ 2 or even smaller), and a sin</?-like cur- 
rent phase relation should be observable. On the other 
hand for large junctions (^4 3> £ 2 ) the statistical fluc- 
tuations become irrelevant, and the behavior should be 
governed by the sin(2<p)-like contribution. As already 
mentioned above a sin(2<^>)-like current-phase relation has 
been observed in experiment 12 . But only some of the 
samples showed this behavior; statistical fluctuations as 
discussed here can be ruled out as a reason, however, be- 
cause the contact in the experiment is quite large com- 
pared to the coherence length of YBCO, £0 ~ 15 A: 
A w 1 /im x 100 nm. A reasonable explanation might be 
a slight deviation from the 45°-orientation of the tilted 
order parameter in some samples, such that the contri- 
bution 1\ is allowed by symmetry. Another explanation 
might be a faceted interface, so that for each facet the 
contribution 1\ is finite. 



IV. CONCLUSION 

To describe a contact between two superconductors 
within the quasiclassical theory, we adapted the bound- 
ary conditions suggested in 2 ^. For numerical calculations 
the choice of a discrete grid of directions is of partic- 
ular importance; we have shown how this grid must be 
chosen in order to guarantee current conservation per- 
pendicular to the interface for an arbitrary scattering 
matrix. Thereafter we proposed a phenomenological ran- 
dom scattering matrix in order to describe microscopi- 
cally rough interfaces without regular structure. 

In section lnTl we applied our model to mirror junctions 
and 45° asymmetric junctions with interface roughness. 
At first we studied the evolution of the temperature- 
dependent critical current for mirror junctions when in- 
creasing the interface roughness: We found that the non- 
monotonic behavior for a specular interface turns to a 
monotonic temperature dependence for rough interfaces. 
Therefore interface roughness is a possible explanation 
for the experimental results. 

For the asymmetric junction we calculated the aver- 
age current-phase relation, which describes large contacts 
We find a leading sin(2(^)-like behavior of the 
current phase relation. This is in agreement with ex- 
perimental findings 1 ^ as discussed in Sec. IIII Bl On the 
other hand we also showed that in small enough junctions 
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a dominating sin ip- like contribution might be present. It 
would be interesting to check this effect experimentally. 

In conclusion our results explain several aspects of the 
experimental data. In order to improve the theoretical 
description, a detailed description of the interface is nec- 
essary in the first place. For example facets, or a modified 
electronic structure near the interface, should be taken 
into account in a more realistic theorj^. 
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